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A scalar-tensor theory of gravity can be made not only to account for the current cosmic ac- 
celeration, but also to satisfy solar-system and laboratory constraints, by introducing a non-linear 
derivative interaction for the scalar field. Such an additional scalar degree of freedom is called 
"Galileon." The basic idea is inspired by the DGP braneworld, but one can construct a ghost-free 
— model that admits a self-accelerating solution. We perform a fully relativistic analysis of linear 

, perturbations in Galileon cosmology. Although the Galileon model can mimic the background evo- 

^ ■ lution of standard ACDM cosmology, the behavior of perturbation is quite different. It is shown that 

JH ' there exists a super-horizon growing mode in the metric and Galileon perturbations at early times, 

suggesting that the background is unstable. A fine-tuning of the initial condition for the Galileon 
fluctuation is thus required in order to promote a desirable evolution of perturbations at early times. 
Assuming the safe initial condition, we then compute the late-time evolution of perturbations and 
discuss observational implications in Galileon cosmology. In particular, we find ant icorrelat ions in 
the cross-correlation of the integrated Sachs- Wolfe effect and large scale structure, similar to the 
normal branch of the DGP model. 
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Q ! I. INTRODUCTION 



^ ■ The origin of the current accelerated expansion of the universe is one of the biggest mysteries in modern 
cosmology. The conventional explanation is that it is caused by a cosmological constant or dark energy, but then 
this idea indicates that the universe today is mostly filled with some unknown energy- momentum component. An 
alternative idea which has been explored actively in recent years is that general relativity (GR) fails to hold at 
cosmological distances, though GR is definitely the most successful theory of gravitation at shorter distances. 

Typical models of modified gravity possess an extra propagating scalar degree of freedom, and the prototype 
, Lagrangian is often described by the Brans-Dicke theory [2|, |3|. However, the theory in its simplest form clearly 
■ contradicts solar-system and laboratory experiments except in the uninteresting parameter range where it is essentially 
indistinguishable from GR In more elaborate models, the extra scalar degree of freedom is carefully designed so 
that the standard A-Cold Dark Matter (ACDM) model is mimicked while solar-system and laboratory tests are safely 
' evaded. In some of /(i?) gravity models Q and their equivalent description in terms of the "chameleon" scalar-tensor 
, theory this is made possible by utilizing a matter-dependent effective potential to hide the scalar field at short 
^ ' distances. Another interesting wav of suppressing the scalar field on short scales is found in the Dvali-Gabadadze- 
. ^ , Porrati (DGP) braneworld model [3], in which a non- linear derivative self-coupling of the scalar is involved Q. Below 
' a certain distance from the source on the DGP brane, the non-linear effect of the brane bending (corresponding to 
I the extra scalar mode) plays an important role even for weak gravitational fields, successfully screening the scalar- 
mediated force • Recently, a general class of scalar-tensor theories has been explored that includes non- linear 
derivative interaction terms in the Lagrangian and shares the same properties as the DGP effective theory [ll| . That 
is, the theory leads to Lorentz invariant equations of motion having only second derivative of the scalar field on a 
flat background. Such a scalar fleld is dubbed a Galileon. The Galileon theory has been covariantized and considered 
in curved backgrounds in Ref. 12] . Cosmology in Galileon modified gravity has been first studied in Ref. [13J. It 
is known that the self-accelerating solution in the DGP model suffers from the presence of a ghost, and therefore it 
is interesting to investigate the ghost issue also in the Galileon-type scalar-tensor theories. In Ref. [l3], a ghost- free 
Galileon model that admits a self-accelerating solution has been constructed. 
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In this paper, the evolution of cosmological perturbations in Galileon modified gravity is studied. We develop a full 
cosmological perturbation theory at linear order, with which we can discuss the stability of the model and compute 
observational signatures such as the integrated Sachs- Wolfe (ISW) effect and weak lensing. These observations help 
us to distinguish between different dark energy/modified gravity models. Though our analysis is not complete in the 
sense that the important non-linear effect is missed, the linear theory is still valid on sufficiently large scales which 
we will be interested in. 

The theory of gravitation we consider is described by the action [ll| - [T^ 



S ■ 



d 



(1) 



where (V(/))^ := g'^'^V^^/'Vi/^, □</> := g^^V v^^: and Cm is the matter Lagrangian. The matter Lagrangian does not 
depend on the Galileon field </>. The most general form of the Galileon action has been derived in Ref. [l^l, but in 
this paper we focus on the simple case which includes only D^CV^)^. Variation with respect to the metric gives the 
Einstein equations 



-U^3m^Va [/(</>) {V4>? 



UJ 

7^ 



(2) 



and variation with respect to the Galileon field leads to the equation of motion 



R + uj 



2a(f> (V0)' 



2/(0) V^V.^V^ V-^/. - (□() 



(V0)' (V0)' + 4^V^(/.V,(/.V^V'^0 = 0. 



(3) 



The model that is most closely related to the DGP braneworld is given by cj = and f{(j)) = r^Mpj/(/)^, where Mpi is 
the Planck scale and is the crossover scale which is taken to be of the order of the present Hubble scale to account 
for the late-time accelerated expansion of the universe [l3| . In this paper we are mainly interested in the case with 
/((^) = rl/cjp and general uj [IJ], but nevertheless the basic perturbation equations are provided without specifying 
the form of /(</>). 

The paper is organized as follows. In the next section, we review the background evolution of Galileon cosmology 
and emphasize the presence of the self-accelerating solution. In Sec. Ill the evolution of linear perturbations in 
Galileon cosmology is addressed analytically and numerically. Observational signatures in the Galileon model are 
discussed in Sec. IV. Finally, we draw our conclusions in Sec. V. 



II. THE BACKGROUND EVOLUTION 



We start with reviewing the background evolution of Galileon cosmology. For the Friedmann-Robertson- Walker 
background metric, ds^ = —dt^ + a'^{t)Sijdx'^dx^ , the Einstein equations give 

3i?^ = ^ - 3HI + + 0V (3ff - f /) /^ (4) 

^3H' -2H = I- + i + I^ + 2HI+^I^-cp^f(i+ /^ (5) 

20 2 \ 2 J 

while the equation of motion for the Galileon field reduces to 
6 (2H^ + - UJ (21 + + 6HI^ 

-6(j)^f {2Hi + HI + 2HI^ + 3H^I^ I + A^'^fail^I + (j)^ f {af + 3ai + aa) = 0, (6) 

where we defined i/(t) — d/a, I{t) := 0/0, and ck„[0(t)] := d"ln//dln0". In the Jordan frame the energy-momentum 
conservation law remains the same as the standard one: p + 3H{p + p) =0. 



3 




FIG. 1: Dimensionless physical distance for different u. 



Let us consider a specific model with 

m - f ' (7) 

for which a self-accelerating solution has been shown to exist At early times, ^ r^, we have an approximate 
solution for Eq. ©: 

Since P ^ <C , the Galileon field is constant on a cosmological time scale, (p ~ I/IGttG, and the Einstein 
equations (HJ and ^ reduce to the usual ones: iH^ — 8nGp and 3H^ + 2H — —SnGp. This is the cosmological 
version of the Vainshtein effect, by which GR is recovered below a certain scale. At late times we find an accelerating 
solution satisfying 



Note that in order for the self-accelerating solution to exist we require 4 -I- 3a; < 0. Interestingly, though 4>is a, ghost in 
this parameter range in the usual Brans-Dicke theory, the Galileon theory does not suffer from ghost-like instabilities. 

It is easy to solve numerically for the intermediate regime between the early-time solution and the late-time ac- 
celerating one. We can confirm that the early-time solution specified above indeed evolves into the self-accelerating 
solution at late times, mimicking the standard ACDM model [14|]. In this ACDM-like solution, the parameter Tc must 
be tuned so as to satisfy, say, Jim («)[:= p/i^'t'H'^)] — 0.3 at the present day. In Fig. [T] we show the dimensionless 
physical distance, r{z) :— Hq dz' /H{z'), for different w as a function of redshift z {:— 1/a— 1). Except for small 
(say, < 10), the background evolution in Galileon cosmology is almost indistinguishable from the ACDM model. 
We shall see how this degeneracy is disentangled at perturbative order. 

III. COSMOLOGICAL PERTURBATIONS 

Let us study the linear evolution of scalar cosmological perturbations. We work in the Newtonian gauge and 
consider a perturbed metric 

ds^ = -(1 + 2$)dt2 + a^(l - 2*)(5jjdxMx^'. (10) 
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The perturbation of the Gahleon field is written as 

d^^^itMt,x), (11) 

while the perturbed energy-momentum tensor is given by 

STo'> = -Sp, ST^^'^d^Sq, ST/^SpS^. (12) 

Here and hereafter an anisotropic stress perturbation is neglected. Since the energy-momentum conservation holds in 
the Jordan frame, we have 

Sp + 3Hidp + 6p)-3'^{p + p) + ^dq = 0, (13) 
Sq + 3HSq+ip + p)<^> + Sp = 0, (14) 
where = S^^didj. The traceless part of the («, j) component of the Einstein equations reduces to 

= (15) 
from which one can see that the Galileon field gives rise to an effective anisotropic stress. It is convenient to define 

f]:=i(vl/ + ci>), (16) 

because weak lensing and the ISW effect are determined through this combination of the metric potentials. Using 
this variable one may write = ^ -\- Lp/2 and $ = Vt — Lp/2. The perturbed Einstein equations and the perturbed 
equation of motion for the Galileon field are summarized in Appendix A without specifying the function /(</>). We 
will, however, focus on the specific model with f{4>) = f'^/4>'^ in the bulk of the present paper. 



A. Early-time behavior 

At early times the effect of the Galileon in the background can be neglected: ^ ^ r~^. Let us look at 
the behavior of cosmological perturbations in this regime. Keeping leading order terms in the Einstein equations, we 
obtain 

6H [n + Hn) -2—n + 8TrGSp = IGnGpip + rll^ [-9Hip+ —if) , (17) 



-2[n + Hn] -SirGSq ^ r'^J^{<f-3Hip), (18) 



TrGSp = -IGnGpip-ririf, (19) 



while the equation of motion for the Galileon becomes a massless wave equation 



0+{3H+^\ip-^^ip^O. (20) 



It is clear from Eq. ([50]) that the fluctuation oscillates inside the sound horizon: 

where we have moved to the Fourier space by promoting to — fc^. Here Co and 9o are the integration constants. 
On super-horizon scales we ignore the spatial derivative term and obtain 

where Ci and C2 are the integration constants. Suppose that the universe is dominated by a single fluid with the 
equation of state w. It follows that the second term grows like ^ ^210/(1+10) (respectively ~ \nt) for it; > (respectively 
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w = 0). This growing behavior is caused by the term {H/H)(p in Eq. (pO)) . which decreases the effect of the Hubble 
friction. The growth of (p then induces the growth of the metric potential on super-horizon scales: 

The growth of the long wavelength modes suggests that the homogeneous background we are considering is unstable. 
The logarithmic growth in the matter-dominated era is rather mild, but in the radiation-dominated stage the C2 mode 
grows like ^ t^/^, which would be more dangerous. Even if C2 = initially so that the above growing solution is 
avoided, the Ci mode influences the evolution of the comoving density perturbation defined by pAc := Sp — 3Hdq. 
Indeed, combining Eqs. ([T7| and ([T8)) one obtains the generalized Poisson equation, 

Ac = {2n + rll^p) ~ 2r2/2-| + (2 - Sr^/^) (24) 

and for long wavelength modes one sees that the last term governs the evolution of the comoving density perturbation: 
Ac ~ constant oc Ci . This will drastically change the growth of the density perturbation, as will be seen in a numerical 
example. 

Since the evolution equation for ip is decoupled from the Einstein equations, we can have an early-time solution 
= by tuning the initial conditions so that ifi = Q and (/j^ = at some t = ti. With this initial conditions, the 
evolution of the metric and density perturbations is identical to that in the ACDM model at early times, and the 
above unwanted behavior could thus be removed. This seems to be the most safe choice of the initial conditions. 
When the effect of the Galileon becomes large in the background and the universe enters the accelerating stage, the 
fluctuation ip will be excited and consequently the evolution of the metric and density perturbations will be modified. 
From now on, we will be focusing on this late-time effect. 



B. Quasi-static approximation 



For the subhorizon evolution of the density perturbation in the matter-dominated stage, the quasi-static approx- 
imation is often employed in the literature, by which the relevant equations can be much simplified. In Galileon 
cosmology, the closed evolution equation for the density perturbation in the quasi-static approximation was already 
derived in Ref. Here, we rederive the equation just for completeness. Later we will confirm the validity of the 
approximation on subhorizon scales by numerical calculations. 

Making the approximation ©(V^^/a^) 0{H^^) inside the Hubble horizon, the conservation equations (fT3)) 
and IHll reduce to 



Ac + 2HAr 



-3>. 



(25) 



(This equation holds exactly on all scales in the ACDM model thanks to the Einstein equations.) Picking up the 
spatial derivative terms in the modified Einstein equations and the Galileon equation of motion, we find 



2r2V^ pAc 

20 



{m -t- rli" (1 + rlf-)\ + 2 (1 + rlf-^ 



-n 



0, 
0, 



(26) 
(27) 



where l3{t) := 3 + 2a; + [41 + 2P 
comoving density perturbation. 



SHI — r^I ) . Using these equations, we obtain the evolution equation for the 



Ac + 2iJ Ac 



pAc 



m 



(28) 



and the modified Poisson equation. 



V2 



m 



pAc 



(29) 
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FIG. 2: Linear growth factor for different k. The initial condition for the numerical calculation is f2i = 1, S7i = 0, and 
Lfi — = ifii. The plotted value is divided by the initial value of the growth factor. The parameters are given by = —50 and 
fimo = 0.3. Results obtained from the quasi-static approximation and in the ACDM model are also shown. 



Since ~ (IBttG)^^ and /3 ^ H/I S> 1 at early times, one can confirm that (ys ~ and the ACDM result is 
recovered. In other words, the quasi-static approximation will no longer be valid in the general early-time evolution 
with iy9 7^ ^ (/?. 

We also have the Galileon fluctuation in terms of the density perturbation: 



P—^^-{l + ril')L^. (30) 



Equation ((30)) would have possible non-linear terms of order rliy"^ / a^Y . Such terms can be neglected provided that 
V'^ip/a? < P/rl- Since Ac (/3/i?^)(V^(/?/a^), this condition is translated to Ac < fi'^ /(jcHY ~ 1- Non-linear terms 
such as r1\/'^ip/a?{yip/aY and r1{yip/a)^ are suppressed by a factor ip{<^ 1). Therefore, the linear approximation is 
valid as long as Ac < 1. The non- linear effect would suppress the deviation from GR via the Vainshtein mechanism. 
Non-linear screening of the scalar degree of freedom is an important aspect of the Galileon-type theories, but it is 
beyond the scope of the present paper. 



C. Numerical solutions 



We are interested in the matter dominant era and the accelerating stage after that, so that we neglect the pressure 
of matter: p = = 5p. Then, the trace part of the Einstein equations and the equation of motion for the Galileon 
imply 



2r2 + {8H + 31- rll^)tl - 2rlPm 



3 + 2w 



lif + ril 



I(p+[2I 



Lp + Hip 



= 



(31) 



and 



(3 -I- 2uj) 



V2 ■ 

(p + iHlp 



6{2H^ + H)-uj(2i + I^ + 6HI] {2n~ip) + 2 



3n + {15H - 



il2H + hl)l(p - [l2Hi + 12ij/ + 27i?/2 + 36i?2/ + 10// + 6/^) ^ -f- ^4/ + 3/^ + 8/// 



-ip 



-2/ 



+ f 6/ -h 10/2 + 2iHi] + I'^n 



-n 



(32) 
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FIG. 3: Evolution of the potential Q for different k. The same parameters as in Fig. [2] are used. 
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FIG. 4: Evolution of the Galileon fluctuation (p for different k. The same parameters as in Fig. [2] are used. 



Given the background evolution and the initial conditions for D, and if, the above two equations can be solved 
numerically in the Fourier space. The matter density perturbation on comoving slices, Ac, is then computed by 
combining the (0, 0) and (0, i) components of the Einstein equations as follows: 



V2 

2—n 



2(f) (f> 



P 



-3/ 



H.I 



2uj 



3m + {%H + AI)m - [6H + -l) ip+ —ip - 



9H^ + -HI + 1^ l^p 



(33) 



This may be regarded as the generalization of the Poisson equation. 

Numerical solutions of Ac, 51, and (p are shown in Figs. [2H7] for uj = —50 and cj = —5000 as functions of redshift z. 
For the present-day value of r2ni(a) we used Unio = 0-3. In order for the perturbation evolution in the early stage to 
be the same as in the ACDM model, the initial condition (pi — — ipi has been used. We find that the quasi-static 
approximation works well for the modes with k > 0.01/iMpc~^, while the evolution deviates from the quasi-static 
solutions on larger scales. As clearly seen, the behavior of the growth rate and metric potential differs from the 
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FIG. 5: Same as Fig. [2] but with lu = -5000. 
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FIG. 6: Same as Fig. El but with u = -5000. 

ACDM case, which will lead to different observational signatures. In particular, it is interesting to note that the 
metric potential J7 grows even when the growth rate decays. The late-time growth of is also found in the normal 
(i.e., non-self-accelerating) branch of the DGP model ^1^. 

To complete the discussion in Sec. IIII we illustrate the effect of the different choices of the initial condition just 
by showing the example obtained from the initial condition i7i = 1, 17^ = 0, Lpi — 0.1 x Vli, and ipi — aX ai ~ 
in Fig. m One sees the oscillating behavior of and tf for short wavelength modes. On large scales the comoving 
density perturbation does not grow, as explained earlier. 

IV. OBSERVATIONAL CONSEQUENCES 

We separate the evolution of the metric and density perturbations from the initial amplitude of $7 and write 
ri(A;; rf) = r])il,{k; rji) and Ac(fc; i]) — 2?A(fc; v)^{k; r]i), where rj is the conformal time, rj = J dt/a, and rji is some 

initial time. The angular power spectrum between two observed fields X and Y on the sky is given in term of the 
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FIG. 7: Same as Fig. g] but with uj = -5000. 
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FIG. 8: Evolution of cosmological perturbations with the initial condition S7i — 1, tpi — 0.1, and tli = tpi = at ai = 10 
The parameters are given by oj = —50 and QmO ~ 0.3. The plotted value is Ac divided by its initial value. 



angular transfer function lf'^{k) by 



-A^(fc)T2(fc) 



if{k)ir{k). 



(34) 



Here, is the primordial spectrum of the comoving curvature perturbation, TZ, and during matter domination we 
have 17 = —(3/5)7?. because the fine-tuned initial condition ipi = Q = ipi implies that the perturbation evolution before 
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the last scattering {rj = rj^,) is identical to the standard one in the ACDM model. The transfer function T{k) is to be 
defined shortly. We consider weak lensing, the ISW effect, and large scale structure (LSS) in the following discussion. 
For weak lensing we have 

lT"^{k)^ r d^-^^Vn{k;rj)Mk{Tjo~v)], (35) 

where ji is the spherical Bessel function, 770 is the present time, and the window function is defined as 

W^(r/) := 2(7?o - f drj' ^l^n{f^'), (36) 
Jo Vo - V 

with n{ri) being a normalized distribution of the source. Similarly, for the ISW effect is given by 

/pW(fc) ^ 2 / drjdMk;v)ji[Hvo - V)]- (37) 

For LSS we have 

rvo 

ll^Hk)^bJ dzng,{z)VA[k;v{z)]Mk{iio-v)], (38) 

where bi is the galaxy bias and ngi{z) is the redshift distribution of the galaxies normalized to J dz ngi{z) — 1 for each 
redshift bin i. In modified gravity, the growth rate is scale-dependent, which would result in the scale-dependence of 
the bias. However, as seen from Figs. [2] and [5l the growth rate is dependent only weakly on scales. Therefore, we 
assume that the bias is scale independent for each redshift bin. 

In the actual numerical calculation we use the conventional scale-invariant primordial spectrum, A^(fc) ~ 2.4x 10~^, 
and the BBKS transfer function [Ig], 

T{k) = ^"^^34^"^ [1 + ^-^^^ + (16-1'^)' + (5-46'?)' + (e-^l?)'] , (39) 

where q = fc/F/iMpc"^ with F = f2,„o/iexp(— f^bo — Vs/i^lbo/^^mo) .17], ^mo = 0.3, fibo = 0.04, and h = 0.7. 

To make the numerical calculations simpler by reducing the number of integrals, we use the approximate formula [isl - 

mi: 

Fix)Ukx)dx ^ ^^IlllI^Fiil + l/2)/fc]. (40) 

The idea that underlies the derivation of the above formula is that the spherical Bessel function ji (x) monotonically 
grows from at a; = to a; ~ Z -|- 1/2 and then starts to oscillate rapidly. For Z ^ 1 we have Z -I- 1/2 ~ I and 
r(Z/2 + 1/2) /Til /2 -I- 1) ~ ^2/Z, and Eq. pO)) results in the Limber approximation often used in the literature. The 
usual Limber approximation is accurate for I > 10, but the approximation at lower / is improved by employing the 
formula pO]) . Indeed, the approximate results obtained from Eq. (|40l) are remarkably accurate even for small Z, as is 
confirmed by comparing them with the exact numerical integration of (k). 
First, we compute the weak lensing convergence power spectrum defined by 

which is shown in Fig. |9l Here we have assumed that the source distribution is given simply by n{r]) — S {r] — ri{zs)) 
with the source location z^ = 3. Although the background evolution is almost the same, the convergence power 
spectrum in Galileon cosmology is enhanced compared to that in the ACDM model, as is expected from enhanced 
n obtained in Sec. IIII CI It is also found that the convergence power spectrum is insensitive to the Brans-Dicke 
parameter uj. As has been done in many other studies [22r.-32j (see, however, Ref. [33l|). the present analysis ignores 
the non-linear growth of structure and non-linear recovery of GR, which would be important for probing modified 
gravity with ongoing and future lensing surveys. The non-linear effect will be more significant at higher multipoles 
(Z> 200). 

Next, let us discuss the ISW effect in Galileon cosmology. The ISW effect is an additional cosmic microwave 
background (CMB) temperature anisotropy on large scales caused by the imbalance between blueshifts experienced 
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FIG. 9: Convergence power spectra for the Galileon models with uu = —20, —50, —5000, and the ACDM model. The source 
location is given hy Zs = 3. 



by the CMB photons falling into the potential well and redshifts imprinted as they climb out. The effect signals a 
net change of the gravitational potential along the path of the photons due to dark energy or the modified gravity 
effect that mimics it. Therefore, it is one of the most useful probes for dark energy and modified gravity (see, e.g., 
Refs. d^-dll for the ISW effect in modified gravity models). We show the ISW power spectra in Galileon cosmology 
in Fig. [TOl The ISW power spectrum is larger in Galileon cosmology than in the AGDM model because of the larger 
time variation of the metric potential Q in Galileon cosmology. It can be seen that the amplitude is larger for smaller 
\Ld\. Note, however, that whether J7 typically grows or decays is not evident from the ISW auto-correlation. (Growing 
decreases the CMB temperature in Galileon cosmology, while decaying increases the temperature in the ACDM 
model.) 

Unfortunately, the ISW effect is typically much smaller than the primary CMB temperature anisotropics and it is 
only significant on the largest scales where cosmic variance is dominant. For this reason, it is difficult in practice to 
extract the ISW effect using CMB observations alone. However, we can measure it by correlating CMB with LSS jiof. 
The ISW-LSS cross-correlation is shown in Fig. 1111 where for definiteness, we assume that the galaxy distribution is 
given by ngi{z) = 1.5 (z^/z'^) exp[— (z/z)^'^] with z = 1. We find anti-correlations in Galileon models as a consequence 
of growing fl at late times, while the cross-correlation in the ACDM is positive due to the decaying potential. The 
Galileon model with smaller |a;| shows a larger anti-correlation signal. Negative ISW-LSS cross-correlations are also 
found in the normal branch DGP model [39]. The cross-correlation between CMB and LSS has been evaluated 
using WMAP data and various different surveys of LSS, and the detection of the positive cross-correlation has been 
suggested [4l| - [5]| . Therefore, the Galileon model is most strongly constrained from the ISW-LSS cross-correlation. 



V. CONCLUSIONS 

In this paper, we have studied cosmological perturbations in Galileon- type modified gravity (TT| - [l3 |. i.e., scalar- 
tensor gravity with a non-linear derivative interaction of the form /((/>) □0(V0)^. The gravitational theory is inspired 
by the effective theory in the DGP braneworld and shares several interesting properties with the DGP model. First 
of all, the Galileon theory admits a self-accelerating solution without ghost-like instabilities, and hence it can be an 
alternative to mysterious dark energy. Moreover, the non-linear derivative interaction suppresses the scalar-mediated 
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FIG. 10: ISW power spectra for the Galileon models with u — —20, —50, —5000, and the ACDM model, obtained by using 
the approximation formula (|40|) . To check the accuracy of the approximation, the result obtained from the exact numerical 
integration of Eq. (|37p is plotted by crosses. WMAP 5-year data for the total CMB angular power spectrum is also plotted 
with error bars. 

force sufficiently near the source, evading solar-system and laboratory tests. This non-linear screening mechanism 
(the Vainshtein effect) works also in a cosmological setting, so that early-time standard cosmology can safely be 
restored. To address the (super-)horizon scale evolution of metric and density fluctuations, we have developed a 
full cosmological perturbation theory in Galileon-type models. It was shown that the perturbation evolution is the 
same as the ACDM case at early times, provided that the initial condition for the Galileon fluctuation ip :— 5(t)/(f) is 
fine-tuned: ipi = Q = (pi. For generic initial conditions we found a super- horizon growing solution in the metric and 
Galileon perturbations, implying that the homogeneous background we are considering is unstable. Even if the initial 
condition is tuned to remove this growing solution, the comoving density perturbation shows an unusual behavior 
on super-horizon scales unless (p — Q ai early times. We thus argue that the most safe initial condition is to choose 
Lpi ^ Q = Lfi aX an early stage in the matter-dominated era. 

We have focused on the evolution of the particular combination of the metric potentials, £7 := + <I>)/2, which 
determines weak lensing and the ISW effect and hence can help to constrain the Galileon model. In the ACDM model, 
f2 (= \[/ = $) decays when A begins to dominate the universe. In contrast to this, it was found that in the Galileon 
model Vt grows rather than decays when the effect of the Galileon field becomes large and the universe begins to 
accelerate. A similar behavior can be found in the normal branch of the DGP model. We have computed explicitly 
the weak lensing convergence power spectrum to show that it is enhanced relative to the ACDM case. The resultant 
power spectrum was insensitive to the Brans-Dicke parameter lo. Note that the effect of the non- linear growth of 
structure must be taken into account in order to predict weak lensing signals correctly. The non-linear effect becomes 
larger at higher multipoles. We have also demonstrated that the cross-correlation of the ISW effect and LSS shows 
an anticorrelation in Galileon cosmology, which will be the most notable signature of the model compared to the 
positive ISW-LSS cross-correlation in the ACDM model. The cross-correlation was negatively larger for smaller |aj|. 
Since positive ISW-LSS cross-correlations have been suggested in the cross-correlation studies of WMAP data with 
different LSS surveys, we conclude that this can provide the most stringent constraint on Galileon cosmology. 

In this paper we have considered the specific form of the non-linear derivative interaction ~ /((/))Dq!)(V0)^ (mainly) 
with /(0) = It would be interesting to explore the effects of higher order interaction terms [Tl|, [T^] on the 

background and perturbation evolution and on the ghost and stability issues. 
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FIG. 11: Cross-correlation between ISW and LSS for the Galileon models with cj = -20, -50, -5000, and the ACDM model. 
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Appendix A: Linearized field equations 

The (00) and (Oi) components of the perturbed Einstein equations are 

6H + H^) - + It - iTT'P = ^^"P - ^HI^ - 3/* - - ujI {ip - /$) 



2(f> 2(f> 



a2 f t2 



fr <j 3/* + 2 {6H - ail) /$ - {9H - 2ail) tp 
V2 



(2 + «,)(3ff-^/U^/ 



-2{^ + H^] = 



Sq 
24> 



2 ) 2 

-{ip - /$) -I- Hp - (1 + Ljj)Iip 

+<l>^fl^ {-/$ + <^ + [(1 + ai)I - 3H] (p} 



(Al) 



(A2) 
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and the trace part of the {ij) equations is 



^' + 3i/* + H^+ i 3H^ + 2i7 $ 



5p p 
20 ^ 20 



/$ + 2 / + r + 2i// $ 



(/3 + 2(i/ + I)(p + ujl{(f- /$) 



//<j -/2$-2 2/+(2 + ai)/2 



-/(^ + 2 



/+ (2 + ai)/2 



(2 + ai)/ + 



2 2 



The traceless part of the (ij) equations simply reads — ^ = ip. The equation of motion for is 



where 



2a; 



<f+{3H + I)ip -ip- (2/ + /2 + 6i?/) $-/($ + 3* 



* + 4i/* + + 4i/^ + 2i/ $ 



(A3) 



V2 

2— (2*-$) 



fi -4:I{3H - ail)<p + 4 -3{HI + HI + SHT + 3H^I + 2ai// + (af + 3ai + 02) / 



ti7(/j + 4 (/ + /2 + 2i// ) —if + 6 



r'J' + 2 // + /^ + 3i/r U' 



-2 (9iJ - 2ai/) + 4 61 [2HI + HI + 2HI^ + 3H^l) - iail^I - {al+3ai+a2) r ^\^0, (A4) 



2/^ — $ 

a"' 

2 f 2 



-6(2 + ai) (^2iJ/ + ij/ + 2i//2 + 3i/2/j / + 4[ai(2 + ai) +02]/^/ 
+ [as + ^2(5 + 3ai) + ai(3 + ai)(2 + ai)] I"^. 
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